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Abstract 

The theory of electron holes is extended into the quantum regime. The Wigner-Poisson system is solved 
perturbatively based in lowest order on a weak, standing electron hole. Quantum corrections are shown to 
lower the potential amplitude and to increase the number of deeply trapped electrons. They, hence, tend 
to bring this extreme non-equilibrium state closer to thermodynamic equilibrium, an effect which can be 
attributed to the tunneling of particles in this mixed state system. 
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Quantum plasmas 1] have recently attracted considerable attention. Non-ideal, dense plasmas 
generated e.g. in the ultraintense laser-solid interaction certainly belong to this category. However, 
also ideal, dilute plasmas -the addressee of this letter- can exhibit a quantum behavior. One reason 
is that the miniaturization of today's micro- and nano-electronic components has reached a level 
such that the system length becomes comparable with the de Broglie wavelength, in which case 
tunneling effects are no longer negligible .2.]. Also, states in combined traps attained by particles 
and their anti-particles used to form anti-hydrogen may be modelled by a plasma having quantum 

n 

features |3|]. Other examples can easily be found, and some of them will be mentioned further 
below. Hence, classical transport models will unlikely be sufficient to describe the plasma behavior 
in such devices adequately. Generally the main focus in this respect is the collective charge carrier 
transport which has been shown jj| to be increased by the presence of coherent structures such 
as hollow phase space vortices j^| . The latter are also called holes in plasma physics Q, 0, 0,13] 
because of their associated density depression exhibiting a remarkable robustness and stability [S ] . 

Before studying the quantum corrections to such holes let us review some further dilute many 
particle systems and related disciplines that reveal quantum aspects. Charged-particle beams 
in particle accelerators are typically dilute systems, so quantum effects are usually disregarded. 
However, a spectrum of phenomena, which recently became more and more important, reveals 
the existence of several quantum aspects of beam physics connecting the physics of particle ac- 
celerator with the frontiers of several disciplines, such as (for instance) plasma physics, radiation 
beam physics, astrophysics, mesoscopic and condensate physics 0, Q|. Most of these phenomena 
introduce a sort of quantum correction to the leading classical behavior of the system. For exam- 
ple, quantum excitation 11] plays a role for the long-term stability of longitudinal electron beam 
dynamics in the high-energy circular accelerating machines while the Sokolev-Ternov effect 
of spin polarization of electron and proton beams is a manifestation, at the macroscopic level, of 
the single quantum nature of the beam particles. Numerical phase space investigations based on 
tracking with the quantum map have shown that quantum corrections can substantially affect the 
particle beam trajectories in the vicinity of the separatrix jl^ ]. 

Recently quantumlike methodologies |lj] have been applied to a number of classical physical 
situations, in which U is replaced by another characteristic parameter of the problem considered. 
For instance, they have been applied to accelerator physics to plasma physics Q|, to surface 
12] and to nonlinear optics 18] in an attempt to describe linear and nonlinear 



gravity wave physics 



problems of the dynamics of beams and large amplitude wavepackets. 

In principle, all these problems can be formulated, in the configuration space, in terms of 
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a system of Zacharov equations, i.e. nonlinear Schrodinger-like equation coupled with one (or 
more) equation(s) taking into account the reaction of the environment. The corresponding phase 
space description is the one provided by the Wigner-Moyal quasidistribution I JJj] whose evolution 
equation, the von Neumann equation, plays the role of a kinetic-like equation associated with the 
system. 

Analytically, the framework under which Wigner-Moyal quasidistributions have been mostly 
considered so far is that of particles interacting with a given external e.g. parabolic potential 
to analyze coherent and squeezed states. Furthermore, a quantum-like phase space analysis of a 
paraxial-charged-particle beam transport, travelling through a quadrupole-like device with small 
sextupole and octupole aberrations, has been carried out showing a satisfactory agreement with 
the results of the standard tracking simulations [^(J and, consequently, the suitability of using the 
quantum phase-space formalism in particle accelerators. This has been done within the framework 
of the thermal wave model Quantum- like corrections involved in the von Neumann equation 
have been discussed for paraxial beams of both particles and radiation and compared with the 
standard classical description j^ll ]. 

Particles in quantum plasmas moving in their own, self-consistent potential, on the other hand, 
have not been given much attention so far. An exception are self-consistent but linearized solutions 
of the Wigner-Poisson system, dealing with quantum corrections to Landau damping of Langmuir 
waves |22| or to the two-stream instability by means of the Nyquist method |23j, and the self- 
consistent linearized solution of the Wigner-Moyal kinetic-like equation for Langmuir wavepackets 



coupled with the ion-acoustic wave equation 



2$ . In particular, the Wigner-Moyal kinetic-like 



description is suitable for describing the Benjamin-Feir instability (modulational instability) as 
well as predicting the stabilizing effect of a sort of Landau damping. It is well known that the 
latter cannot be shown in configuration space, where the the system is usually described by the 
Zakharov equations. By using the pure state formalism, a Landau-type damping has been shown 
for the longitudinal dynamics of both charged-particle coasting beams and e.m. wavetrains in 



high-energy circular accelerators and nonlinear media, respectively 



25j] . A similar approach has 



been extended (mixed state formalism) to ensembles of partially-incoherent waves in different 



physical situations 



24 



26, 



27]. 



The aim of this paper is to describe an electron-ion unmagnetized plasma, in which, on the 
basis of the experimental evidences as well as on theoretical and numerical investigations mentioned 
above, the quantum nature of the particles is not disregarded. However, it is taken into account 
only as a weak (perturbative) effect in comparison to the leading classical behavior of the system. 
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Together with the weak quantum effect, we take into account the usual classical electrostatic 
collective plasma effects coming from the standard meanfield approximation of the coulombian 
interaction, in such a way that our system is described by a set of coupled equation comprising the 
von Neumann equation for the Wigner-Moyal quasidistribution and Poisson's equation. Hereafter, 
we will refer to this system of equations as WP-system (Wigner-Poisson system) . Our goal is to find 
a self consistent solution of the WP-system to the lowest order of the quantum correction. Before 
formulating our problem in detail, it is worth mentioning some further considerations starting from 
the classical case. 

Classically, electron and ion holes are nonlinear, stationary solutions of the Vlasov-Poisson 

mey are 



system (VP-system) being omnipresent structures in many driven, collisionless plasmas. 

29], in particle accelerators jsol. I3I | - in the 



found in one or other variant in the laboratory |2S|, 
laser-plasma interaction and in the extraterrestrial space jfl j ] . Analytical solutions, applicable 
to kinetic hole structures found in particle accelerators, have been presented in 

Generally speaking, the fundamental role of holes arises from the fact that they can nonlinearly 
destabilize a plasma even in linearly stable situations, namely if they posses a negative energy 



character |35 



361 ] . There is hence accrued interest to extend the studies of holes into the quantum 



domain, which means that the WP-system has to be employed instead. 

In recent numerical studies a multistream model for a current-driven quantum plasma 
has been applied. Signatures of coherent hole structures appear in the simulation of a statistical 
mixture of many pure states, with each wave function obeying the Schrddinger-Poisson system. 



Sal. To the best 



That this system is equivalent to the mixed state WP-system has been shown in 
of our knowledge, an analytic self-consistent nonlinear solution of the WP-system is still missing 
in the literature. 

In this letter we present a first rigorous nonlinear self-consistent solution of the WP-system 
assuming weak nonlinearity and proximity to the classical VP-system. First we shall refer to the 
classical e-hole and then study quantum corrections, for which a self-consistent solution is derived. 

We are first investigating a standing, classical e-hole which is the simplest inhomogeneous 
stationary solution of the VP-system for a plasma with immobile ions: 

[vd x + <f/(x)a v ]f(x,v)=0, (la) 



dvf(x,v)-l:=n e (x)-l, (lb) 
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where f(x,v) is the distribution function of electrons in phase space. Here space, velocity and 
density are normalized by the electron Debye length Ad = (fcr e /47rnoe 2 ) 1 / 2 , the electron thermal 
velocity vth = {kT e /m e ) 1 /' 2 and the homogeneous electron density, no, where m e is the mass of the 
electrons and T e is their temperature. 

In thermal equilibrium, the plasma adopts a homogeneous state with a Maxwellian distribution 
in velocities /m(x, v) = -J= exp (— v 2 /2) . Referring to the potential method jsjj, we solve (jlaj) by 



f(x, v) = J= [eM-E)9(E) + eM-PE)9(-E)} , (2) 

V Z7T 



where E = v 2 /2 — (ft represents the single electron energy. The separatrix E = in phase space 
distinguishes the free (E 1 > 0) from the trapped (E < 0) electron population. Note that / in (J5J), 



reduces to Jm as </> — > 0. 
half power expansion of ( 



'he electron density in the weak amplitude regime can be written as a 

n e {<ft) = l + 4>-^ 7 ^4? /2 + 0{ ( ft 2 ). (3) 
Defining the pseudo-potential as —V'(4>) := n e ((ft) — 1, we find from © with F(0) = 



2 V 15-v/i 

It has to fulfill two necessary conditions: (a) V{<ft) <OinO<0< , and (b) = 0, where 

ip is the amplitude of the perturbation in the potential, which is assumed to be small, i/i < 1. 
From (b) we arrive at —(5 = 15 y/n/ 16 y/^jj — 1 ~ 15-y/7r/ 16^/0 ^ 1- Therefore, j3 has to be a large 
negative number, corresponding to a depletion of the distribution in the trapped particle range. 
On the other hand, equation Q allows us to integrate Poisson's equation (jlb|) . and we obtain the 
bell-shaped electrostatic potential: 

(ft(x) = -0sech 4 (^j . (5) 
Note that other electrostatic structures such as propagating electron holes, ion holes or periodic 



nonlinear waves (cnoidal waves) can also be found by appropriate extensions of this method 

y 



To study quantum corrections, we start with Wigner's quasidistribution which satisfies the time 
independent von Neumann or quantum Liouville equation 19]: 



vd x f + — 

l£ 



x+^d v )-cf>(x-^d v )\j 0. (0) 



where e is the dimensionless Planck's constant e = — — 



Assuming that e is small, i.e. the quantum effects appear only as corrections to the classical 
solution, we can perform a power expansion of the potential operator <\> (x±yc\,) which we insert 
into (|H|). All even terms cancel out and we get up to the third order 

vd x f + 4>'{x)d v f - ^'"( x )d 3 v f = 0, (7) 

which is the equation we have to couple with Ijlbj) . 

As we keep terms up to 0(e 2 ), we will look for corrections of the same order in the potential 
and in the distribution function, / = /o + e 2 /i, = </>o + £ 2< Ai> with /q and 4>q representing now 
(J2J) and © respectively. Inserting this ansatz into (J7| and ()lb|) neglecting again terms of 0(e 4 ) we 
find 

[vd x + <t>' (x)d v ] h = -<t>[(x)d v f + -Lfoixy'dtfo, (8a) 



d 2 x Mx)= I dvh. (8b) 

By defining g(x,v) := f\ + (pidEfo, we can reduce equations (|5aj) and (|8b|) to the somewhat 
simpler system 

[vd x + <f>' (x)d v ] g = ±-^{x)dlf Q =: h(v, x), (9a) 



/+oo 
dvg(x,v), (9b) 
-OO 

2 

Now it is convenient to switch into a new set of variables defined by £ = x, E = ^ — </>o (x), 
a = sg(v) and rewrite h(x,v) = H(£, E,a), g(x,v) = G{^,E,a). With these variables, (|9*a|) 
becomes d^G(^, E, a) = H(^, E,a)/v(£,, E, a), whose general solution is 

G&E,a)=G(0,E,a)+ t^ ^f^ (10) 

Jo w (4 ,E,a) 

where v(^,E,a) = cr-v/2 [E + 0q(£)]. Therefore, in order to find G we only have to integrate 
H(^, E, cr)/v(^, E, a) along the classical particle trajectories given by E = const.. In this expression 
we have chosen the lower integration limit as £ = because this is the only point which is reached 
by all trajectories (see below). Note that a trapped particle will move along a closed, bounded 
trajectory around the origin in phase space. 
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Now we need to replace H (£, E, a) by its full expression. Denoting / in Eq. © as fo(E) we get 
by differentiation 

d 3 J = -Li> { [3 - 2(£ + O )] e- £ 0(i?) + (3 2 [3 - 2/3(£ + O )] e"^0(-£) 
V 2vr 

- [3(1 - /3) - 2<^o(l - /? 2 )] - 2(E + to)(l " /?)*'(£)} 
1 (11) 



'2tt 

For positive energies, we can follow the trajectories up to any £ in (|10j) and, assuming that the 
correction vanishes at £ — ► ±oo, we arrive at 

°<°- * ') - /_! < l^r - sts /I * * ({,)fi({ ' £) - (12) 

Note that this expression does no longer depend on er. For negative energies G(0, E, a) is not 
determined by such a procedure but, due to the symmetry of the problem we can assume that it 
will also be cr-independent. On the other side, we can always extend the integration of (|10j) to 
— oo for negative energies also as long as we change the integration constant. Therefore we have, 
for any E, 

1 



G(Z,E,a)=G(£,E) :-- 



G (E) + [ d?<l%gM?,E) 



(13) 



3!4V27T 
with G (E) = for E > 0. 

Note that 8e/o is discontinuous at E = 0. Therefore, G(£,E) does not have a definite value at 
the separatrix. Our approach will be to solve © for positive and negative energies separately and 
then put both solutions together imposing the continuity of f\ at the separatrix. 
In order to integrate (|13|) . we consider these two different cases: 

1. For E > we have G (E) = and Sl(£,E) = [3 - 2{E + MO)} e~ E - The integral (O can 
be performed analytically to yield 

G(£,E) = -±=[M0 2 + ^-2E-2MO)^o(0]e- E . 



2. If E < we must take G${E) into account. In this case Q(£,E) = 
(3 2 [3 - 2{E + MO)0\ and {BJ) reads 

G(£, E) = -±j= {g (E) + [PMO 2 + (3 - 20E - 2(3MO) 4(0] P 2 e~ 0E } ■ (14) 
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FIG. 1: a) Correction to the potential for ip = 0.1; b) the classical and the quantum corrected potential. 

The continuity of f± is now imposed to determine Go(E). As f\ = g — (PiOe/o, the discontinuity 
of g(x,v) = G(£,E), namely AG := G(£,0+) - G(£,(T) should be equal to <t>iA(d E fo). Since it 
holds A(d E f ) = {13- l)/v^F, we get 

AG = siim K(O 2 + (3-20o(£M'(O 

-0 1 [M(e) 2 + ( 3 - 2/300(0) #,'(£)] " G (0)} . (15) 

Then we can find MO as MO = AG/A(d E fo) = V^AG/(p - 1) with AG given by (US). 
Moreover, as we impose 0i(±oo) = 0, we know that Go(0) = 0. Hence, we obtain 4>i(0- Figure 
^ shows 0i(£) for ^ = 0.1 and Fig. represents the corrected potential. For reference, also the 
unperturbed potential is drawn. We see that the potential experiences a reduction as a result of 
quantum correction. 

To determine Go (E) for all negative energies we go back to ()9b|) . As we have already determined 







a V 
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(f>i(x), the left hand side is now given. It is convenient to write it in terms of <f>Q. To do this we 

/ \l/2 

note that all derivatives of cpo can by expressed by 4>o itself via </>q(£) = —<f>o(£) ( 1 — \/ 4>o((,) / ) i 
= O (£) (l - 5vVo(0/W4) and also V"(0o) = - (l - lSy/fa/i/)/^. Inserting these ex- 
pressions into (|9bl) we find an expression for its left hand side as a function of </>o which we call 
L((fto). An explicit form of this function is too long to be included in this letter and will be presented 
elsewhere (4J|. 

On the other hand, the right hand side of l)9b|) can be written as 

+oo ^ ( ^)=j> r dE ^i-=2r dE^i (i6) 

„ J -MO v (^ E > a ) J -MO V K^ E A) 
Making use of (|TT?|) we can reduce lfTfi|) to 

+°° \ q (E) 



Where R((j>o) represents a known function that can be obtained analytically 



To perform the remaining integral of (|17|). we make a half power expansion of Gq(E): 

-±—G (E) = a l/2 \E\ 1 / 2 + ai \E\ + a 3/2 |£| 3 / 2 + . . . (18) 

OlZV Z7T 



With this ansatz, we have 

And we can finally reduce ()9b|) to 



F^r(i + f) 



L(0 O ) = R(<f>0) + V^\/^ E W l i!n! fl n/ 2 ( 2 0) 



2 ^ r (l + t) 

Therefore, if we define p(i) := j (L(t 2 ) — R(t 2 )), we can find all a n / 2 as 

i /2"r(| + f) d>(t) 



a n/2 



n! V vr r (l + § ) dt n 



(21) 



t=o 



With this expression for a n / 2 we can sum Gq(E) and then find G(£,E) and /i. The correction 
of the distribution function, fi(x,v) is plotted in Fig. El while the final, corrected distribution 
function / = /o + £ 2 fi is represented at fixed x in Fig. 03 

We clearly recognize a partial filling of the phase space within the separatrix being maximum at 
the hole center. An interpretation may be given in terms of refraction or tunneling: in the classical 
solution nearby its separatrix, the region of untrapped electrons is populated stronger than that of 
trapped electrons. In the quantum domain when tunneling becomes effective this gives rise to a net 



FIG. 2: Correction of the distribution function fi(x, v). 



influx of particles resulting in a less dilute distribution of trapped electrons. This interpretation 
of the collective particle behavior, found analytically,matches well with the numerical findings of 
Ref. IrI according to which the quantum corrections affect the particle trajectories in the vicinity 



of the separatrix (see also Ref. 



20J). 



We, therefore, conclude that the overall effect of a quantum correction to a classical e-hole is 
the tendency of the system to reduce the coherent excitation by both a diminution of the amplitude 
and a partial filling of the trapped particle region by refraction (tunneling), bringing the system 
closer to the thermal state. 

Open questions are how these semiclassical corrections are modified in case of finite amplitudes 
\E r > 0(1), of finite quantum corrections e > 0(1), of hole propagation vq > 0, of nonlocality 
of structures such as periodic wave trains, some of which will be addressed in our forthcoming 



publication 



4l|. 



This work was performed under the auspices of the DAAD (Deutscher Akademischer Aus- 
tauschdienst) and CRUI (Conferenza dei Rettori delle Universita Italiane) in the framework of the 
research program "VIGONI." 
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